Bound states and scattering lengths of three two-component 
particles with zero-range interactions under one-dimensional 

confinement 

0. I. Kartavtsev and A. V. MalykhQ 

Joint Institute for Nuclear Research, Dubna, 141980, Russia 

S. A. Sofianos 

Physics Department, University of South Africa, Pretoria 0003, South Africa 

(Dated: October 24, 2008) 

Abstract 

The universal three-body dynamics in ultra-cold binary gases confined to one-dimensional motion 
are studied. The three-body binding energies and the (2 + l)-scattering lengths are calculated for 
two identical particles of mass m and a different one of mass mi, which interactions is described 
in the low-energy limit by zero-range potentials. The critical values of the mass ratio m/mi, at 
which the three-body states arise and the (2 + l)-scattering length equals zero, are determined 
both for zero and infinite interaction strength Ai of the identical particles. A number of exact 
results are enlisted and asymptotic dependences both for mjm\ — * oo and Ai —* — oo are derived. 
Combining the numerical and analytical results, a schematic diagram showing the number of the 
three-body bound states and the sign of the (2 + l)-scattering length in the plane of the mass ratio 
and interaction-strength ratio is deduced. The results provide a description of the homogeneous 
and mixed phases of atoms and molecules in dilute binary quantum gases. 
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I. INTRODUCTION 



Dynamics of few particles confined in low dimensions is of interest in connection with 
numerous investigations ranging from atoms in ultra-cold gases to nonos- 



tructures 



lOj] . Experiments with ultra-cold gases in the one-dimensional (ID) and 



quasi- ID traps have been recently performed 



HQ 



12 



est to the investigation of mixtures of ultra-cold gases 



13 1, amid the rapidly growing inter 



17. 



19, 



20|. Different 



aspects of the three-body dynamics in ID have been analyzed in a number of recent papers, 



e. g., the bound-state spectrum of two-component compound in 2l|], low-energy three-body 
recombination in 22], application of the integral equations in {23], and variants of the hy- 



perradial expansion in 



24 



25 



26]. 



It is necessary to emphasize that the exact solutions are known for an arbitrary number of 
identical particles in ID with contact interactions 



27 



281 ]; in particular, it was found that the 
ground-state energy E N of N attractive particles scales as E N /E N=2 = N(N 2 — l)/6. There 



is a vast literature, in which the exact solution is used to analyze different properties of 



ew- 



and many-body systems; few examples of this approach can be found in Ref. |29l . 1301 . 1311 . [32]. 

The main parameters characterizing the multi-component ultracold gases, i. e., the masses 
and interaction strengths can be easily tuned within wide ranges in the modern experi- 
ments, which handle with different compounds of ultracold atoms and adjust the two-body 
scattering lengths to an arbitrary values by using the Feshbach-resonance and confinement- 



resonance technique |33j. Under properly chosen scales, all the properties of the system 
depend on the two dimensionless parameters, viz., mass ratio and interaction strength ratio, 
the most important characteristics being the bound-state energies and the (2 + l)-scattering 
lengths. In particular, knowledge of these characteristics is essential for description of the 
concentration dependence and phase transitions in dilute two-component mixtures of ultra- 
cold gases. 

In the present paper, the two-component three-body system consisting of a particle of 
mass mi and two identical particles of mass m interacting via contact (5-function) inter- 
particle potential is studied. In the low-energy limit, the contact potential is a good approx- 
imation for any short-range interaction and its usage provides a universal, i. e., independent 



of the potential form, description of the dynamics 



23, 



26 



35, 



36. 



37j. More specifically, 



it is assumed that one particle interacts with the other two via an attractive contact interac- 



2 



tion of strength A < while the sign of the interaction strength Ai for the identical particles 
is arbitrary. This choice of the parameters is conditioned by an intention to consider a 
sufficiently rich three-body dynamics since the three-body bound states exist only if A < 0. 

Most of the numerical and analytical results can be obtained by solving a system of 
hyper-radial equations (HREs) 38] . It is of importance that all the terms in HREs are 
derived analytically; the method of derivation and the analytical expressions are similar 



36 



37[. To 



to those obtained for a number of problems with zero-range interactions [26, 
describe the dependence on the mass ratio and interaction-strength ratio for the three- 
body binding energies and the (2 + l)-scattering length, the two limiting cases Ai = 
and Ai — > oo are considered and the precise critical values of m/m\ for which the three- 
body bound states arise and the (2 + l)-scattering length becomes zero are determined. 
Combining the numerical calculations, exact analytical results, qualitative considerations, 
and deduced asymptotic dependencies, one produces a schematic "phase" diagram, which 
shows the number of the three-body bound states and a sign of the (2 + l)-scattering lengths 
in the plane of the parameters m/mi and Ai/|A|. This sign is important in studying the 
stability of mixtures containing both atoms and two-atomic molecules. 

The paper is organized in the following way. In Sect. [IT] the problem is formulated, 
the relevant notations are introduced, and the method of "surface" function is described; 
the analytical solutions, numerical results and asymptotic dependencies are presented and 
discussed in Sect. IIII( the conclusions are summarized in Sect. IIVI 



II. GENERAL OUTLINE AND METHOD 



The Hamiltonian of three particles confined in ID, interacting through the pairwise con- 
tact potentials with strengths Aj, reads 

where X, and are the coordinate and mass of the ith particle, Xjk = Xj — x^, and {ijk} is 
a permutation of {123}. In order to study the aforementioned two-component three-body 
systems, one assumes that particle 1 interacts with two identical particles 2 and 3 through 
attractive potentials and denotes for simplicity m 2 = m 3 = m and A 2 = A3 = A < 0. The 
corresponding solutions are classified by their parity and are symmetrical or antisymmetrical 



3 



under the permutation of identical particles, depending on whether these particles are bosons 
or fermions. The even (odd) parity solutions will be denoted by P = (P = 1). 

In the following, the dependence of the three-body bound state energies and the (2 + 
l)-scattering lengths on two dimensionless parameters m/mi and Ai/|A| will be investigated. 
Hereafter, one lets h = \\\ = m = 1 and thus m\ 2 /h 2 and h 2 /(m\\\) are the units of energy 
and length. Furthermore, one denotes by A and A\ the scattering lengths for the collision 
of the third particle off the bound pair of different and identical particles, respectively. 
The scattering length is considered at the lowest two-body threshold, which corresponds to 



determination of A if Ai/|A| > — y 2/(1 + m/rrii) and A\ otherwise. With the chosen units, 
£?th — — 1/[2(1 + m/mi)] and E' th = — A?/4 are two-body thresholds, i.e., the bound-state 
energies of two different and two identical particles, respectively. 

The binding energy and the scattering length are monotonic functions of the interaction's 
strength and for this reason much attention is paid to calculations for two limiting cases 
of zero (Ai = 0) and infinite (Ai — > oo) interaction between the identical bosons. It is of 
interest to recall here that due to one-to-one correspondence of the solutions [39( all the 
results derived for systems, in which the identical particles are bosons and Ai — > oo, are 
applicable to those in which the identical particles are fermions and the s-wave interaction 
between them is zero (Ai = 0) by definition. 

The —a. and analytical results will be obtained mostly by giving a system of 



HREs 



381 ] where the various terms are derived analytically 26j, |36|, [37|. The HREs are 



written by using the center-of-mass coordinates p and a, which are expressed via the scaled 
Jacobi variables as psina = X2 — £3 and pcosa = cotcj (2x\ — X2 — £3) given the kinematic- 
rotation angle uj = arctan JY 



expanded as in papers 



24 



25 



2m/m-[ so that E, 



26, 



th 



cos 2 co\ The total wave function is 



P 



-1/2 



£/n(p)$n(«>P) 



(2) 



n=l 



in a set of functions $ n (a, p) satisfying the equation at fixed p 

/ d 2 



da'- 



(3) 



complemented by the condition 

d$ n (a,p) 
da 



a=w+0 



a=w— 



2pcosu§ n (ui, p) = 



(4) 



which represents the contact interaction between different particles 2y, |35|, |37|, |40| . Taking 
into account the symmetry requirements, one can consider the variable a within the range 
< a < n/ 2 and impose the boundary conditions 



(1-T)^ + T$ n 
oa 







- a=Tr/2 



o=0 







(5) 
(6) 

oo). These 



where P = (P = 1) for even (odd) parity and T = (T = 1) for Ai = (A x 
boundary conditions are posed if two identical particles are bosons, however, the case T = 1 
is equally applicable if two identical particles are noninteracting (Ai = 0) fermions. 

The solution to Eq. ([3]) satisfying the boundary conditions ([5]) and (J6j) can be written as 

cos[£ n (cj — 7r/2) — Pn/2)] cos(£ n a — Tn/2) , a < uj 
cos(£„c<j — Tn/2) cos[£ n (a — tt/2) — Pn/2] , a > uj 
where the normalization constant is given by 



$„(a,p) = B. n 



(7) 



Bf 



dp 



(8) 



— 2 cos 2 (£ n {c<j — 7r/2} — Ptt/2) cos 2 (£„u; — Tir/2) cosu 

In order to meet the condition (TJJ, the eigenvalues £ n (p) should satisfy the equation 

2p cos uj cos[£ n w - (£„ + P)tt/2] cos(^ - Tvr/2) + £ n sin[(£ n + P- T)n/2] = . (9) 

Notice that the case P = 1 and T = is formally equivalent to the case P = and T = 1 
under the substitution of for tt/2 — uj. 

The expansion of the total wave function (T2]) leads to an infinite set of coupled HREs for 
the radial functions f n {p) 

d 2 C(P) ~ 1/4 



dp 2 



+ E 



oo 

/n(p)"E 
m=l 

nn 



Using the method described in [26|, I36L 
terms in Eq. (flOl) . 



d d 

(P) - Qmn(p) j 3-Qmn(p) 

dp dp 



Up) = o . (10) 

37j, one can derive analytical expressions for all the 



Qnm(p) = ($„ I &J 



PUP) = « I *'J 



Q, 



"n ' m 



+ 







e" W 


5 




~ J 



1 (e' r 



6 £' 4 V £' 



n ^ m 



n = m 



(11) 



(12) 



where e n = ^ and the prime indicates derivative with respect to p. 

The obvious boundary conditions for the HREs ffTUl) f n (p) — > as p — > and p — > oo was 
used for the solution of the eigenvalue problem. For the calculation of the scattering length 
A, one should impose the asymptotic boundary condition for the first-channel function 

~ psina;- A, (13) 

while all other boundary conditions remain the same as for the eigenvalue problem. The con- 
dition (fT3|) follows from asymptotic form of the threshold-energy wave function at p — > oo, 
which tends to a product of the two-body bound-state wave function and the function de- 
scribing the relative motion of the third particle and the bound pair. The linear dependence 
of the latter function at large distance between the third particle and the bound pair leads 
to asymptotic expression ffl3|) for the first-channel function in the expansion (T5]). On the 
other hand, the expression (TlBl is consistent with the asymptotic solution of the first-channel 
equation in fflUj) . in which the long-range terms Pn{p) and — l/(4p 2 ) cancel each other at 
large p. 

III. RESULTS 

A. Exact solutions 

There are several examples, where the analytical solution of the Schrodinger equation 
for the systems under consideration can be obtained. Firstly, for a system containing one 
heavy and two light particles (in the limit m/m\ — > 0), using the separation of variables, 
the solutions can be straightforwardly written both for zero and infinite interaction strength 
between the light particles. In particular, for Ai = 0, there is a single bound state with 
binding energy E% — — 1 and the (unnormalized) wave function is 

^ b = e -I^H*i3l ; ( 14 ) 
whereas the scattering wave function at threshold energy E t ^ = —1/2 is 

^sc = (M - 1) e-'* 131 + (|x 13 | - 1) e-l^l , (15) 

which gives the (2 + l)-scattering length A = 1. On the other hand, for Ai — > oo, the 
three-body system is not bound, and the scattering wave function at the threshold-energy 

6 



E th = -1/2 is 

V sc = \x 12 e-^-x 13 e-^ \ , (16) 

which gives A = 0. 

Furthermore, as mentioned in the introduction, the exact solution is known for an 



arbitrary number iV of identical particles with a contact interactions in ID [27J, [28] 
and if the interaction is attractive there is a single bound state, which energy equals 
= —N(N 2 — l)/24. In particular, for three identical particles (m = m\ and Ai = A) 
there is only one bound state with energy E% — — 1 and the (unnormalized) wave function 



is 



*b = exp(--X;|xy|l , (17) 
whereas the exact scattering wave function at the two-body threshold E t ^ = E' th = —1/4 is 

^ sc = ^exp(-^|x 4i |) -4exp(-i^|a; ii |) , (18) 

i<j i<j 

which implies that the (2 + l)-scattering length is infinite \A\ — > oo, i. e., there is a virtual 
state at the two-body threshold [24]. 

Further exact results can be obtained by using the abovementioned correspondence of the 
three-body solutions for the infinite interaction strength (Ai — > oo) between two identical 
bosons and for two noninteracting fermions (Ai — > 0). For example, for three equal- mass 
particles (m — mi) the exact wave function at the two-body threshold (E& = —1/4) reads 

f e -xis/2 + x 12 /2 _ 2e-*23/2 x > o 

* SC= , ,2 /2, (19) 

^ | e «Dls/2_ e xi 2 /2| j Sis < . 

As follows from (|T9|) . the (2 + l)-scattering length is infinite; as a matter of fact, this 



implies a rigorous proof of the conjecture 



211 ] that m = mi is the exact critical value for the 



emergence of the three-body bound state in the case of infinite repulsion (Ai — > oo) between 
two identical bosons. 

It is worthwhile to recall here the exact solution for three equal- mass particles (m = m±) 
if the interaction between two of them is turned off (Ai = 0) 4l(. A transcendental equation 
was derived for the ground-state energy, which approximate solution gives the ratio of three- 
body and two-body energies E 3 /E th f& 2.08754. 



7 




FIG. 1: Mass-ratio dependences for the even-parity states; shown are the ratio of the three-body 
bound-state energies to the two-body threshold energy (left) and the (2 + l)-scattering length A 
(right). Presented are the calculations for a system containing two identical bosons with zero (solid 
lines) and infinite (dash-dotted lines) interaction strength Ai. The dash-dotted lines represent also 
the results for a system containing two identical noninteracting (Ai = 0) fermions. Encircled are 
those points, in which the exact analytical solution is known. 

B. Numerical calculations 

For the even-parity states (P = 0) and the two limiting values of the interaction strength 
between identical bosons, Ai = and Ai — > oo, the HREs (jlOp are solved to determine the 
mass-ratio dependence of three-body binding energies and the (2 + l)-scattering length A. 
The calculations show sufficiently fast convergence with increasing the number of channels; 
15-channel results are presented in Fig. [TJ The precise critical values of the mass ratio, for 
which the three-body bound states arise (\A\ — > oo) and the (2 + l)-scattering length A = 
are presented in Table d and are marked by crosses in Fig. [1] and Fig. [3j The condition that 
the ground state energy is twice the threshold energy is important as it determines whether 
production of the triatomic molecules is possible in a gas of diatomic molecules. The mass 
ratio, at which E^/E t h = 2 is determined to be m/mi ~ 49.8335 for Ai — ► oo, while for 
the excited states the condition Es/E t h = 2 is satisfied for m/m\ ~ 130.4516 if Ai = and 
mjm\ ~ 266.1805 if X 1 — ► oo. 
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TABLE I: The even-parity critical values of the mass ratio m/m\ for which the (2 + l)-scattering 
length becomes zero (marked by A = 0) and an nth three-body bound state arises (marked by 
| A | — > oo). Calculations done for two values of the interaction strength between the identical 
particles, Ai = and Ai — > oo. 







Ai 


= 




Ai 


— ► oo 


n 


m/mi{A = 


0) 


77i/?7li(|j4| — ► Oo) 


m j yn\ (A = 


0) 


7TI / 77t^(L4 — > Oo) 


1 








0* 




1* 


2 


0.971 




2.86954 


5.2107 




7.3791 


3 


9.365 




11.9510 


16.1197 




19.0289 


4 


22.951 




26.218 


32.298 




35.879 


5 


41.762 




45.673 


53.709 




57.923 


6 


65.791 




70.317 


80.339 




85.159 


7 


95.032 




100.151 


112.179 




117.583 


8 


129.477 




135.170 


149.222 




155.193 


9 


169.120 




175.374 


191.463 




197.989 


10 


213.964 




220.765 


238.904 




245.973 



Exact 



As shown in Fig. [lj the binding energies increase with increasing the mass ratio, whereas, 
the scattering length A has a general trend to decrease with increasing the mass ratio on each 
interval between two consecutive critical mass ratios at which the bound states appear. Nev- 
ertheless, the calculations for Ai = show that A{mjm\) becomes non-monotonic function 
at small m/m\. More precisely, the scattering length takes a maximum value A « 1.124 
at m/mi Pd 0.246. Again one has to note that the mass-ratio dependence of energy and 
scattering length (plotted in Fig. [1]) and the critical values of the mass ratio (presented in 
Table [I]) are the same both for the three-body system containing two identical bosons if 
\\ — ► oo and for the three-body system containing two identical noninteracting (Ai = 0) 
fermions. 

It is of interest to note that the calculated binding energy Es/E t h ~ 2.087719 for three 
equal-mass particles (m = mi) if two identical ones do not interact with each other (Ai = 0) 
is very close to the result 41] E 3 /E t ^ w 2.08754 obtained from the analytical transcendental 
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FIG. 2: Mass-ratio dependence of the (2 + l)-scattering length A for odd-parity states (P = 1) of 
a system containing two identical noninteracting bosons (Ai = 0). The numerical calculation (solid 
lines) is compared with the large-mass-ratio asymptotic behaviour given by Eq. (|25p (dash-dotted 
lines). The dependence corresponding to large A > 15 is shown on a large scale in the inset. 



equat 



ion (see Sect. IIII A|) . A small discrepancy most probably stems from the approximations 



of |4l| made in numerical solution of the transcendental equation. The (2 + l)-scattering 
length turns out to be small and negative, A m —0.09567, for m = mi and Ai = and takes 
a zero value at slightly smaller mass ratio m/m\ ~ 0.971 (see Tabled]). 

Analogously, the odd-parity (P = 1) solutions for three-body system containing two iden- 
tical noninteracting bosons (Ai = 0) were obtained. As follows from Eq. (Q, the eigenvalues 
£n(p) entering in HREs ({TO], are nonnegative, which implies that there are no three-body 
bound states. The calculated dependence of the scattering length A is shown in Fig. [2l 
A increases monotonically with increasing mass ratio following the asymptotic dependence 
discussed in Sect IIII CI 



C. Asymptotic dependencies 

1. Large attractive interaction of two identical particles 

In the limit of large attractive interaction between the identical particles, Ai — > — oo, the 
even-parity wave function takes, with a good accuracy, the factorized form ~ <t>o{x23)u(y) 
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[y = cotu;(2xi — £2 —£3)], where (fio(x) = J\xjj2 exp(— |Aix|/2) is the wave function of 
the tightly bound pair of identical particles with energy E' th = — A 2 /4 and u(y) describes the 
relative motion of a different particle 1 with respect to this pair. Within this approximation, 
u(y) is a solution of the equation 



d 2 

— + 2|Ai| exp 1 



1 + 2m/m 1 \Xiy\) + A 2 /4 + E 



u(y) = 



(20) 



which gives the independent of Ai leading-order terms in the asymptotic expansion for the 
three-body binding energy, e ~ 4/(1 + 2m/mi), and the (2 + l)-scattering length, 



Ai « (1 + 2m/ mi) /A ■ 



(21) 



2. Two heavy and one light particles 

For large mass ratio m/rrix, one can use the adiabatic and quasi-classical approximations 
which provide, e. g., a universal description for the energy spectrum j^j. To describe the 



three-body properties in the limit of largem/m! — > 00 [to — > tt/2— Jmx/(2m)\, one considers 
the first eigenvalue £i(p) = in(p), which large-p asymptotic dependence is approximately 
given by 

pcosu; = ; 77; 7 — r-^- , (22) 

H 1 _|_ (_l)P e -K(7r-2w) ' ^ ) 

as follows from Eq. on the equal footing for the system containing two identical bosons 
both for Ai = and Ai = 00 and for the system containing two identical noninteracting 
fermions. 

The number of the three-body even-parity (P = 0) bound states n can be determined 
for large m/mi, using the one-channel approximation in ffTUj) and the effective poten- 
tial —K 2 (p)/p 2 , from (I22p . Within the framework of the quasi-classical approximations 
and taking into account the large-p asymptotic dependence (|22|) . one obtains the relation 
m/mi ~ C(n + 5) 2 in the limit of large n and mjvn\. The constant C can be found as 

C=— / V2t + t 2 ) -r^O 

2 [Jo 2t{l + t) 2 

where the integral is expressed by letting t = exp[— k(7t — 2a;)] in the leading term of the 

quasi-classical estimate, 



2 

w 2.59 , (23) 



COS LO I dp 
JO 
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2 



1/2 



1} = Tin . (24) 



Fitting the calculated mass-ratio dependence of the critical values, at which the bound states 
appear, to the n- dependence C(n + 5) 2 (up to n — 20, see Table fl] for 10 lowest values), 
one obtains in a good agreement with the quasi-classical estimate (l2"5|) C ~ 2.60 both for 
Ai — ► oo and Ai = 0. Simultaneously, one obtains 5 = 0.73 if X\ — > oo and 5 = 0.22 if Ai = 
for the parameter, which determine the next-to-leading order term of the large-n expansion. 

The asymptotic dependence of the effective potential —n 2 (p)/p 2 obtained from (1221) allows 
one to find the leading order mass-ratio dependence of the odd-parity (P = 1) scattering 
length, 



where 7 ~ 0.5772 is the Euler constant. The convergence of the calculated dependence 
A{mjm\) to the asymptotic dependence (1251) is shown in Fig. [2]for the case of two identical 
noninteracting bosons (Ai = 0). 

D. Mass-ratio and interaction-strength ratio dependencies 

Collecting the numerical and the exact analytical results, the asymptotic expressions, and 
qualitative arguments, one obtains a schematic "phase" diagram, which depicts the number 
of three-body bound states and the sign of the (2 + l)-scattering lengths in the m/m\ - 
Ai/|A| plane (shown in Fig. [3]). 

The plane of parameters is divided into two parts by a dotted line, Ai/|A| = 



— y2/(l + m/mi), with the low-energy three-body properties being essentially different in 
the upper and lower part, where the two-body threshold is determined by the bound-state 
energy of two different and identical particles, respectively. The lines which represent the 
condition \A\ = 00 or \Ai \ = 00 (arising of the three-body bound state) separate areas with 
different number of the bound states, whereas the conditions A = or A\ = split each 
area into two parts of different signs of the scattering lengths. 

It can be proven rigorously that in the upper part of the diagram (above the dotted line), 
the number of the three-body bound states n increases and the (2 + l)-scattering length 
A decreases with decreasing the interaction strength Ai, while in the lower part (below the 
dotted line) n increases and Ai decreases with decreasing the mass ratio m/mi. The proof 
is based on the representation for which the lowest two-body threshold is independent of 
Ai and mi in the former and latter case, respectively. The required conclusion follows from 




(25) 
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/ \ 1/2 
(m/mj 



FIG. 3: Schematic "phase" diagram for the even-parity states of two identical bosons and the 
third different particle. The dotted line marks the border between two areas where the lowest 
two-body threshold is set by the energy of two different and two identical particles. The number 
of the three-body bound states is marked by n in the corresponding areas separated by solid lines. 
The sign of the (2 + l)-scattering lengths A and Ai is marked by db and the corresponding areas 
are separated by dashed lines. The crosses show the calculated critical values of the mass ratio 
(enlisted in Tabled]). Encircled are those points, in which the exact analytical solution is known. 

the monotonic dependence of the Hamiltonian on Ai and mi. A schematic "phase" diagram 
demonstrated in Fig. [31 is drawn by using more strict assumption on the positive slope of 
the lines, which show where the three-body bound states arise (\A\ — > oo) and where the (2 
+ l)-scattering lengths (A = and A\ = in the upper and lower part of the Ai/|A| - m/rax 
plane, respectively) become zero. Tentatively, this assumptions seems to reflect correctly 
the general trend; nevertheless, one should note that the slope of the isolines of constant 
scattering length is not generally positive. In particular, A is not a monotonic function of 
the mass ratio for Ai = 0, as shown in Fig. [2j this implies a non-monotonic dependence of 
the constant-A isolines in a region near the point (m/mi = 0, Ai/|A| = 0). 

For sufficiently large repulsion Ai and small mass ratio m/mi the three-body bound 
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states are lacking. The limit m/m\ — > (ID analogue of the helium atom with contact 



interactions between particles) was discussed in paper 42j, where the binding energy as a 
function of the repulsion strength between light particles was calculated and the critical value 
of the repulsion strength for which the three particles becomes unbound was determined. 
Recently, a very precise critical value Ai/|A| ~ 2.66735 was found in [2l|. The boundary of 
the n = area (shown in the upper left corner in Fig. [3]) goes from the point (m/mi = 0, 
Ai / 1 A | « 2.66735) to the point (m/mi = 1, Ai — > oo), as was conjectured in [21] and proven 
in Sect. IIII AI by using the exact solution at the latter point. Taking into account this 
result, the above- discussed monotonic dependence on Ai, and the exact solution for three 
identical particles, one comes to an interesting conclusion that there is exactly one bound 
state (n = 1) of three equal-mass particles independently of the interaction strength Ai. 
There is exactly one bound state (n = 1) also for a sufficiently large attraction between 
identical particles whereas the second bound state appears for m > m x and |Ai| < 1 (as 
shown in Fig. [3]). Therefore, the scattering length A± changes from the positive value given 
by (12 II) at Ai — > — oo to the negative one as Ai increases. The strip areas corresponding to 
n > 1 are located at higher values of the mass ratio with the large-n asymptotic dependence 



n oc ym/rrii. In each parameter area corresponding to n bound states, the scattering lengths 
run all the real values tending to infinity at the boundary with the n — 1 area and to minus 
infinity at the boundary with the n + 1 area. 

IV. CONCLUSION 

The three-body dynamics of ultra-cold binary gases confined to one- dimensional motion is 
studied. In the low-energy limit, the description is universal, i. e., independent of the details 
of the short-range two-body interactions, which can be taken as a sum of contact 5-function 
potentials. Thus, the three-body energies and the (2 + l)-scattering lengths are expressed 
as universal functions of two parameters, the mass ratio m/mi and the interaction-strength 
ratio Ai/|A|. The mass-ratio dependences of the binding energies and the scattering length 
are numerically calculated for even and odd parity and the accurate critical values of the mass 
ratio, for which the bound states arise and the scattering length became zero, are determined. 
It is rigorously proven that m/mi = 1 is the exact boundary, above which at least one 



bound state exists (as conjectured by [211]) : the related conclusion is the existence of exactly 
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one bound state for three equal-mass particles independently of the interaction strength 
between the identical particles. Asymptotic dependences of the bound-state number and the 
scattering length A in the limit m/m\ — > oo and of the binding energy and the scattering 
length Ai in the limit Ai — > — oo are determined. Combining the numerical calculations, 
analytical results, and qualitative considerations, a schematic diagram is drawn, which shows 
the number of the three-body bound states and the sign of the (2 + l)-scattering length as 
a function of the mass ratio and interaction-strength ratio. 

The obtained qualitative and quantitative results on the three-body properties provide 
a firm base for description of the equation of state and phase separation in dilute binary 
mixtures of ultra-cold gases. In particular, a sign of the (2 + l)-scattering lengths essentially 
controls the transition between the homogeneous and mixed phases of atoms and diatomic 
molecules. The condition E^/E^ > 2 defines the parameter area, where the production of 
the triatomic molecules is energetically favorable in a gas of diatomic molecules. 

From the analysis of the "phase" diagram in Fig. [3] it follows that still there are inter- 
esting problems deserving further elucidation. These include the problem of non-monotone 
dependence of the constant-A isolines in the Ai/|A] - m/rrii plane, the behaviour of the 
lines separating the positive and negative scattering lengths within the n = 1 area, and the 
description of the beak formed by the lines separating the n = 1 and n = 2 areas in the 
vicinity of the exact solution for three identical particles (Ai = A and m = mi). 

One should discuss the connection of the present results with those, which take into 
account the finite interaction radius R e and (quasi)- ID geometry. The determination of the 
corrections due to finite interaction radius is not a trivial task, however, one expects that 
the corrections should be small for all calculated values provided R e /a and R e /a\ are small, 
where a and a\ are the two-body scattering lengths. On the other hand for sufficiently tight 
transverse confinement, one expects that the main ingredient is the relation between the 
3D and quasi-lD two-body scattering lengths established in 33J]. Moreover, a role of the 
transverse confinement does not simply reduce to renormalization of the scattering lengths; 
the full scale three-body calculations are needed to determine the energy spectrum and the 
scattering data in the (quasi)- ID geometry. 

It is worthwhile to mention that more few-body problems are of interest in binary mix- 
tures. In particular, the low-energy three-body recombination plays an important role in 
the kinetic processes, while the elastic and inelastic cross sections for collisions either of 
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diatomic molecules or of atoms off triatomic molecules are needed to describe the properties 
of the molecular compounds. 
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